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INTRODUCTION  AND  SUMMARY 


Aerojet  ElectroSystems  Company  has  been  performing  research  for  the 
Office  of  Naval  Research  on  a  number  of  topics  in  electromagnetic  theory  of 
interest  in  explaining  the  microwave  radiometric  characteristics  of  sea  ice. 
Under  Contract  N00014-81-C-0553  (see  Report  6140  published  in  December  1981 
and  the  Final  Report  7788  published  in  March  1985)  the  bilocal  approximation 
in  strong  fluctuation  theory  was  developed  and  applied  to  isotropic  random 
media  such  as  snow.  The  work  reported  on  here  and  funded  under  Contract 
N00014-83-C-0726  is  an  extension  of  these  ideas  to  the  anisotropic  case 
which  is  characteristic  of  most  forms  of  sea  ice.  Specificall^  this  study 
has  concentrated  on  deriving  equations  for  computing  the  tensor  dielectric 
constant  of  sea  ice  and  making  comparisons  of  the  computations  with  published 
measurements  on  various  components  of  the  dielectric  tensor  (as  well  as 
attenuation  coefficients,  which  are  closely  related  to  the  dielectric  constant) 
The  theory  is  shown  to  be  able  to  account  quantitatively  for  larqe 

parts  of  the  published  experimental  data  on  a  variety  of  ice  types  in  the 

.C 

frequency  range  0.1  to  40  GHz  and  temperatures  from  -5  to  -306C.  A  few  dis¬ 
crepancies,  however,  are  noted.  It  is  not  known  whether  further  development 
of  the  theory  will  be  required  to  resolve  these  differences  between  theory 
and  experiment  or  whether  the  experiments  themselves  must  be  refined.  It  is 
well  known  that  accurate  dielectric  constant  measurements  on  sea  ice  present 
some  rather  difficult  experimental  problems.  Some  inconsistencies  within  the 


existing  body  of  experimental  data  itself  are  pointed  out  in  this  work 
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AM  ANALYSIS  OF  THE  TENSOR  DIELECTRIC  CONSTANT  OF  SEA  ICE 
AT  MICROWAVE  FREQUENCIES 

I.  INTRODUCTION 

Sea  ice  is  an  interesting  and  complex  substance  of  considerable 
importance  for  remote  sensing  studies  of  the  earth.  An  essential  ingredient 
of  any  attempt  to  achieve  a  quantitative  understanding  of  the  microwave  radio- 
metric  and  radar  signatures  of  the  diverse  types  of  naturally  occurring  sea 
ice  is  a  knowledge  of  the  dielectric  properties  of  the  ice.  A  summary  of 
experimental  work  on  the  dielectric  constant  involving  both  artificial  and 
natural  sea  ice  may  be  found  in  [1].  To  this  may  be  added  the  work  of 
Bogorodskii  and  Khokhlov  [2]  which  emphasizes  work  performed  in  the  Soviet 
Union.  Further  experiments  are  currently  underway.  However,  the  highest 
frequency  which  appears  to  have  been  investioated  experimentally  is  still 
below  the  40  GHz  limit  already  discussed  in  [1J  and  [2]. 

The  attempts  to  analyze  the  results  of  these  experiments  have 
ranged  from  linear  regression  models  [13,  [3],  Wiener's  dielectric  mixina 
formula  [3]  and  a  simple  averaging  formula  for  spherical  inclusion  [4],  to 
more  sophisticated  models  which  explicitly  recognize  the  fact  that  sea  ice, 
at  least  in  some  of  its  forms,  should  be  anisotropic  due  to  the  elonoated 
brine  cells  which  it  contains.  In  the  latter  cateaory  are  the  formulas 
based  on  Vant's  dissertation  discussed  in  [1]  and  a  variant  due  to 
Farrelly  [5]. 

All  of  these  formulas  are  strictly  limited  because  scatterina 
effects  due  to  the  random  brine  cells  and  air  pockets  in  the  ice  are  not 
treated.  Thus,  for  example,  they  are  not  well  suited  for  studyino 
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multi-year  ice  at  frequencies  above  20  GHz  because  of  the  dominant  role 
played  by  scattering  from  air  bubbles-  It  is  the  intent  of  this  work  to 
extend  the  methods  of  strong  fluctuation  theory,  which  has  already  been 
developed  for  isotropic  random  media  [6],  to  calculate  the  dielectric  con¬ 
stant  of  an  anisotropic  medium  such  as  sea  ice.  By  using  this  approach, 
scattering  effects  are  included.  Improvements  compared  to  previously  calcu¬ 
lated  behavior  are  obtained  even  at  frequencies  as  low  as  0.1  GHz.  This  is 
partly  due  to  the  fact  that  the  strong  fluctuation  result  for  the  quasi¬ 
static  dielectric  constant  does  not  agree  with  Tinga's  formula  which  was  used 
by  Vant  and  Farrelly  and  partly  the  result  of  the  greater  flexibility  which 
stronq  fluctuation  theory  allows  in  specifying  the  statistical  properties  of 
the  brine  cells.  It  is  not  limited  to  assuming  elliptically  shaped  cells  as 
was  necessary  in  [1]  and  [5]. 

II.  PHYSICAL  MODEL  OP  SEA  ICE 

The  model  to  be  discussed  here  assumes  that  sea  ice  is  a  mixture  of 
pure  ice,  brine,  and  air.  The  pure  ice  occurs  in  the  form  of  thin  platelets 
of  the  order  of  0.5  mm  thick  with  the  crystalloqraphic  c-axis  perpendicular 
to  the  plane  of  the  platelets  [7].  These  platelets  are  stacked  together  to 
form  grains  whose  average  cross  sectional  diameter  is  about  0.4  cm  near  the 
surface  of  the  ice  and  increases  linearly  by  1  cm  for  every  30  cm  of  depth 
[8).  The  c-axis  of  the  grains  can  have  various  orientations.  For  frazil, 
the  axis  may  be  assumed  to  be  randomly  oriented  |1]  while  for  columnar  ice 
there  is  a  tendency  for  the  polar  anqle  of  the  c-axis  (measured  from  the 
vertical)  to  increase  from  approximately  0°  near  the  surface  to  90°  as  the 
depth  increases.  The  azimuthal  orientation  of  the  c-axis  vas  taken  to  be 
random  in  this  study.  This  imposes  a  certain  symmetry  on  the  dielectric 
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tensor  which  will  hold  in  most,  but  not  all,  cases  of  interest  and  excludes 
ice  formed  in  the  presence  of  strongly  directional  currents. 

The  brine  cells  in  the  model  were  assumed  to  occur  between  the 
platelets  comprising  a  grain  and  are  thus  perpendicular  to  the  c-axis. 
Structures  such  as  brine  drainage  channels  were  ignored.  Within  a  grain  the 
brine  cells  occur  as  long,  roughly  cylindrical,  parallel  inclusions  with 
random  spacings.  The  radius  of  these  cells  is  of  the  order  of  0.025  mm  [7] 
while  their  length  varies  from  3  mm  [1]  to  over  1  cm  [7], 

For  purposes  of  the  model,  air  pockets  in  the  ice  are  assumed  to 
occur  only  between  the  ice  grains.  Thus,  the  grains  consist  of  only  brine 
and  ice.  For  lack  of  more  definitive  information  the  air  bubbles  are  assumed 


to  be  roughly  spherical  in  shape  so  that  they  may  be  described  by  a  spheri¬ 
cally  symmetric  correlation  function. 

In  order  to  calculate  the  dielectric  properties  of  sea  ice,  it  is 
necessary  to  specify  the  dielectric  constants  of  the  constituents  as  well  as 
their  respective  volume  fractions.  For  this  study  the  dielectric  constant 
equations  discussed  in  [9]  for  pure  ice  and  in  flO]  for  brine  in  sea  ice 
were  used  while  the  dielectric  constant  of  air  was  taken  to  be  one.  Assuninc 
a  mean  density  Pice  =  0.926  g/cm^  [8J  for  bubble-free  sea  ice,  the  volume 
fraction  of  pure  ice  is 

vice  =  ^  P  Pbvb'//Pice 

where  p  is  the  density  of  the  sea  ice  sample  under  consideration,  the 

density  of  brine,  and  vp  the  volume  fraction  of  brine.  The  volume  fraction 
of  air  is 
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vair  “  1  ~  vb  _  vice  ( 2 ) 

A  fit  to  the  data  of  Nelson  [11]  yields 

pb  «  1-02814  -  0.88128xl0"2t  -  0.9298xl0~4t2  9/cm3  (3) 

in  the  temeprature  range  -1  >t  i  -54°C.  The  brine  volume  may  be  expressed 
as 

vb  =  Sp  (4) 

where  S  is  the  ice  salinity  (absolute,  not  parts  per  thousand)  and  p  is 
a  function  of  temperature  which  has  been  tabulated  in  [8]-  For  numerical 
work  the  fit 

-2.28  -  52-56  t_1  if  t  »  -2.06°C 

0.930  -  45-917  t_1  if  -2.06  >  t  >  -8.2°C 

1.189  -  43.795  t_1  if  -8.2  >  t  >-22.9°C 

21.9921  +  2. 968 56x10 3 t-1  +  1 . 53039xl05t-2  +  3. 502798xl06t-3 
I  +  3.0401xl07t“4  if  ~22*9  >  t  >  -36.8°C 

2.8167  +  0.9494xl0-1t  +  0.9603xl0_3t2 

if  -36.8  >  t  >  -43 . 2°C  (5) 

is  convenient.  The  above  equations  for  p  in  the  range  t  >  -22.9°C  were 
first  published  by  Frankenstein  and  Garner  (12). 

III.  BILOCAL  APPROXIMATION  FOR  THE  DIELECTRIC  CONSTANT 

The  use  of  the  bilocal  approximation  in  strono  fluctuation  theory 
for  a  random  medium  yields  the  equation  [13] 


Report  7975 


Vx  V  x  Em  -  kg  V0Em  -  ■  0  (6) 

for  the  mean  electric  field  propagating  in  the  medium.  Here  k0  is  the  free 
space  propagation  constant,  KQ  the  quasi-static  dielectric  tensor  and  £(r) 
a  random  tensor  defined  by 

Mr)  =  AM  I  +  S  AL]-1  (7) 

where 

AL  =  kg(K0  -  Kr  (r) )  (8) 

The  tensor  S  is  the  coefficient  of  the  delta  function  part  of  the  Green's 
function  F  which  satisfies  the  equation 

V  x  7  x  r  -  kg  Kor  =  I  6 (r  -  r‘ )  (9) 

and  I  is  the  unit  dyad.  Kr(r)  in  (8)  is  the  random  dielectric  tensor  of 

the  medium.  The  last  term  in  (6)  accounts  for  the  effects  of  scattering  from 
the  random  inhomogeneities.  Its  i*-'1  component  is  given  by  the  principal 
value  integral 

t(|G|)Em)i  =  PV  [dV  (4ij(r)  tkl<£'  ))  r}k(r,r'  )  ( r '  )  (10) 

where  F'  is  the  kernel  of  G  and  is  related  to  F  by 

F*(r,r')  =  F( r , r ' )  -  S  6(r  -  r ' )  (11) 

In  (10)  and  below,  the  summation  convention  is  used  for  repeated  Latin  letter 
subscripts  and  the  annular  brackets  denote  expected  value. 

If,  as  in  [6],  a  plane  wave 

Ep  =  E0  p  exp( i  k  •  r  ]  (12) 
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with  polarization  p  is  considered  in  a  homogeneous  medium  and  boundary 
effects  are  neglected  so  that  r  becomes  an  infinite  medium  Green’s  function 
whose  argument  depends  only  on  the  difference  r-r*  ,  it  is  found  that 
satisfies  the  wave  equation 

V  x  V  x  Ej  -  kj  K  E”  =  0  (13) 

where  the  effective  dielectric  constant  K  is 


o 


The  scattering  contribution  F  in  (14)  is  defined  by 

Ei£  =  pvjd3r’  Cijk^(r  -  r' )  rjk(r  -  r’ )  exp[ -ik • ( r-r ' ) )  (15) 

where 

Cijkf(£-E’>  =  <4ij<£>  4k £<£’>)  <16> 

A  translation  of  the  origin  shows  that  E  is  actually  independent  of  r  . 

Specific  equations  for  the  components  of  KQ  and  E  will  be 
developed  later.  However,  on  the  basis  of  the  physical  model  discussed 
above,  the  general  structure  of  the  dielectric  tensor  K  can  be  described. 

If  the  z  axis  is  chosen  to  be  vertical,  the  random  azimuthal  orientation  of 
the  c  axes  of  the  ice  grains  implies  that  K0  is  invariant  under  rotations 
about  the  unit  vector  z  and  reflections  in  any  plane  containing  the  z  axis. 
Thus  Kc  is  necessarily  a  diagonal  tensor  of  the  form 

Kc  =  diag(  K01 ,K0i ,K03 )  ( 1 7  I 

The  situation  for  K  is  somewhat  more  complicated  because  of  the  occurrence 
of  the  propagation  vector  k  in  (15).  If  k  were  zero,  then  K  would  have 
the  same  structure  as  K0  .  However,  since  we  are  interested  in  finite 
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particle  size  to  wavelength  ratios,  additional  components  are  not  excluded  in 
K  -  namely,  additive  terms  proportional  to  kk  ,  kz  ,  and  zk  each  of  which 
are  form  invariant  under  the  above-mentioned  transformations.  One  may  con¬ 
clude  that  the  coefficients  of  kz  and  zk  are  identical  because  E  in  (15) 
is  a  symmetric  tensor  as  a  result  of  the  symmetry  of  C  in  its  first  two  and 

last  two  indices  and  the  fact  that  the  infinite  medium  Green's  function  is 

A  A 

symmetric.  Let  the  y  axis  be  chosen  so  that  k  lies  m  the  x-z  plane: 

k  =  (kx,0,kz )  (18) 

Then  these  considerations  show  that  K  is  of  the  form 


K  *  I  0 


0  Ki 


Koo  0 


k13  0  k33 


where  each  of  the  components  are  functions  of  kx  and  ky  .  Furthermore, 

K^3  vanishes  if  either  kx  or  kz  in  (18)  is  equal  to  zero. 

In  tyoical  applications,  one  of  the  components  of  k  is  known. 

The  value  of  K(k)  may  be  computed  after  solving  a  dispersion  equation  for 
the  unknown  component  of  k  .  The  latter  is  obtained  from  substitution  (12) 
into  (13).  If  p  =  y  ,  then 

k2-k^!<22  ( k )  =  0  (20) 

The  solution  of  this  equation  yields  the  effective  dielectric  constant  for 

,  .  'N  A  A 

the  ordinary  wave.  If  p  lies  in  the  x-z  plane,  the  extraordinary  wave  is 
determined  by  the  determinanta 1  equation 


r  2  .  1,2  _ 


-  K] 


-k  k  -  k2K,  , 
x  z  o  1  3 


-k  k  -  k2K  , 
x  z  o  1 3 


k2  -  k2  -  k 2  K 

z  o  33 


r*  -•  .*  ■ 

•  a  *•  *  m  4 
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To  carry  out  this  program,  it  remains  to  specify  the  components  of 
K  more  explicitly.  The  quasi-static  and  scattering  contributions  to  K 
for  sea  ice  are  treated  below. 

IV.  THE  QUASI -STATIC  DIELECTRIC  CONSTANT 

In  the  model  of  sea  ice  the  quasi-static  dielectric  constant  K0 
is  computed  in  two  stages.  First,  an  expression  is  found  for  the  equivalent 
dielectric  constant  of  the  grains  described  in  Section  II  and  secondly  this 
expression  is  used  to  compute  the  quasi-static  dielectric  constant  of  the  ice 
as  a  whole  taking  into  account  the  arbitrary  orientations  of  the  c  axes  of 
the  grains  and  the  air  bubbles. 

For  an  individual  grain,  let  us  choose  a  coordinate  system  x'  , 
y'  ,  z'  oriented  so  that  the  brine  cells  are  parallel  to  z'  .  The  quasi¬ 
static  dielectric  tensor 

K'  =  diag(K1,K1,K3)  (22) 

of  this  grain  is  determined  according  to  strona  fluctuation  theory  by  the 
equation  [6] 

^  J grain  ^  =  ®  (23) 

where  tqrain  is  determined  by  equations  of  the  form  (7)  and  (8)  with  K0 
replaced  by  K'  and,  since  the  grain  is  assumed  to  be  a  two-component  medium 
consisting  of  pure  ice  and  brine,  the  random  tensor  Kr  reduces  to  a  r  ndom 
scalar  multiple  of  the  unit  dyad  which  may  assume  the  values  K^ce  or 
Kbrine  (the  dielectric  constants  of  ice  and  brine  respectively).  Thus,  in 
(23),  the  expected  value  reduces  to  a  sum  of  two  terms  weiqhted  according  to 
vice'/  *  vb+vice '  and  vb^vb+vice*  which  are  the  respective  probabilities  of 
a  point  being  in  ice  or  brine. 
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The  tensor  Sg  to  be  used  in  (23)  and  corresponding  to  S  in  (7) 
depends  on  the  geometry  of  the  brine  cells.  Specific  formulas  are  discussed 
in  114].  Because  of  the  axial  symmetry  of  the  brine  cells,  Sg  is  diagonal 
of  the  form 


Sq  =  diagtSgj,  Sg2,  Sg3) 


(  24) 


with 


sg3  =  -<*o2  +  2K1Sgl)/K3 


(25  ) 


If  the  equi-correlation  surface  is  assumed  to  be  defined  by  the  ellipsoid 


(x'2+y'2)/£2  +  z'2/#2  =  constant 


(  26  ) 


then 


*gi 


■[(1  +  e2)tan_1e  -  e]/(2k£l<3e3) 


(2-7) 


where 


e  =  [ (^2K3)/(£2K1)  -  1 ] 1/2 


(28) 


This  choice  corresponds  somewhat  to  the  model  of  elliptic  inclusions  devel¬ 
oped  in  [1]  and  [5]  if  the  correlation  lengths  Ip  ,  and  (?z  are  assumed  to 
be  proportional  to  the  radius  and  lenqth  of  a  mean  brine  cell.  It  will  be 
seen  later  that  this  is  not  the  most  satisfactory  choice  for  the  equi- 
correlation  surface.  In  fact,  the  roughly  cylindrical  shape  of  the  brine 
cells  suggest  that  an  equi-correlation  surface 

(x '  2+y '  2  )  l/2/(?p+  |z  '(  /£z  =  constant  (29) 

may  be  more  realistic.  Using  the  formulas  of  [14],  this  assumption  yields 

sgi 


-  j^-2a-1/2  £n[  a*/2  -  1]  -  1  +  [a  -  1]1/2J^’Z/ 


[2 lp  a  k2(K1K3)l'/2] 


( 30  ) 


CL 
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where 

a  =  1  +  (K1C|)/(K3e2)  (31) 

The  study  in  [15]  indicates  that  iz  and  £p  should  be  proportional  to  the 

length  and  radius  of  a  mean  brine  cell.  On  the  basis  of  the  data  in 

Section  II,  the  ratio  £z/£p  should  be  larger  than  100. 

After  K'  has  been  determined  according  to  the  above  procedure,  a 
% 

second  application  of  the  prescriptions  of  strong  fluctuation  theory  leads  to 

the  quasi-static  dielectric  constant  of  sea  ice  which  is  treated  as  consisting 

of  ice  grains  and  air  bubbles.  To  accomplish  this,  it  is  necessary  to  use 

the  standard  x,y,z  coordinate  system  oriented  with  the  z  axis  vertical. 

In  this  coordinate  system,  the  grain  dielectric  constant  becomes  K£  which 

ft, 

is  related  to  £'  by 

k£  =  A-1  K'  A  (32) 

~9  R5  2=  ftj 

where  A  is  the  matrix  transforming  from  the  standard  coordinate  system  to 
ft: 

the  prime  system.  Introducina  Euler  anqles  6B  ,  4>B  ,  *  where  the  conven¬ 

tions  of  [16]  are  used,  and  notinq  that  yB  may  be  set  equal  to  zero  because 
of  the  assumed  axial  symmetry  of  the  brine  cells  about  the  z'  axis,  A  is 
given  by 


cos<|>B 


sin<t>B 


A  =  I  -cos0Bsin4>B  cos0Bcosd>B  sin0B 
\  sin0Bsin<i>B  -sin0Bcos4>B  cos0B 


where  6b  is  the  angle  between  the  axis  of  the  brine  cells  in  a  grain  and 
the  vertical  and  4>b  describes  the  azimuthal  position.  Note  that  £5  is 
random  because  these  angles  are  random. 


in 
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The  equation  ^  =  0  which  determines  the  quasi-static  dielectric 

ft 

constant  [6]  of  sea  ice  may  be  written 


=  0 


(34) 


vair  £air+vqrain  4g 

where  va^r  is  given  by  (12)  and  Vgrajn  is  the  sum  of  and  Vice  . 

These  quantities  are  the  probabilities  of  a  point  being  in  the  air  bubbles  or 
grains  respectively.  The  tensors  £  in  (34)  are  of  the  form  shown  in  (7) 
and  (8).  In  the  case  of  ,  Kr  =  I  in  (8).  Here,  and  below,  the  symbol 

—  represents  an  average  over  the  angles  0B  and  4>B  so  that  |q  is  the 

<«w  “ 

angular  average  of  £a  which  is  determined  by  (7)  and  (8)  if  (32)  is  used 
for  Kr  .  The  air  bubbles  are  assumed  to  be  roughly  spherical  and  hence 
described  by  a  spherically  symmetric  correlation  function.  Thus  S  will  be 
a  diagonal  tensor 


£  =  diag ( Sj , S3 , S3 ) 


(35) 


of  the  form  (25)— (28)  with  the  substitutions  K3  —  Kqi  ,  K3  — *Kq3  ,  and  the 
restriction  H p  =  Pz  .  With  this  expression  for  S  a  routine  computation 

shows  that  the  components  of  tQ  are 

- 

4 ql  1  =  j  °1  Y  -  *3*3  4  (  S1S3Pcos24>B  +  OsSiCOS2^  -  0383) 

Psin20B|/(S1o1V)  (36a) 


4q22  =  |°lY  -  °'i°3  +  ( S3S3P  sin24>B  4  <»3S3Sin24>B  -  O3S3] 

psin26B  |/(  SiO’jY ) 


( 36b ' 


^33  =  !y  -  +  S1psin2eB|/(S3V) 


(36c) 


11 
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£gl2  =  4g21  =  Pla3  +  s3 PI  sin20B  cos4>B  sin4>B/(ci'1,y )  (36d) 

4gl3  *  Ig3i  *  ~Pcos6B  sin0B  sin<t>B/Y  (36e) 

Ig23  =  4g32  =  pcos6B  sin0B  cos4>B/y  (36f) 

where 

Oi  =  1  +  k2Si(Koi  -  K*)  li  =  1  or  3)  (37a) 

P  =  k2(K3  -  K2)  (37b) 

y  =  +  (O1S3  -  -  S3S3 p)psin20B  (37c) 


Part  of  the  average  to  obtain  £g  fron  these  equations  is  easy  to  perform. 
The  assumption  of  equally  probable  azimuthal  orientations  of  the  c  axis  is 
equivalent  to  integrating  over  $B  and  dividing  by  2~  .  This  immediately 
yields  a  diagonal  tensor 

4  g  =  ^ia9 ( 4gl * 4gl > 4g3 ^  (37) 

However,  the  0B  average  depends  on  introducing  a  probability  distribution 

P (0 B )  for  the  brine  axis  angle  before  integrating.  For  all  types  of  ice 

except  frazil,  which  will  be  treated  below,  a  lack  of  knowledqe  of  P(0B) 

precludes  a  rigorous  integration.  To  circumvent  this  difficulty  the  mean 

value  theorem  may  be  applied.  This  is  equivalent  to  treating  0B  as  a 

parameter  representing  an  appropriate  average  brine  cell  angle.  Thus  £n 

ft 

may  be  evaluated  from  (36a)-(36c)  after  replacing  cos2<5B  and  sin2(fcB  by 

1/2. 


12 
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For  frazil,  where  all  orientations  of  the  c  axis  are  equally 

probable,  P(0)  =  1/2.  In  this  case,  the  complete  symmetry  requires  Kn  ,  S  , 

5SJ  SS 

and  to  be  scalar  multiples  of  the  unit  dyad  (denoted  by  K0f  ,  Sf  and 

*** 

4f  respectively.  Remembering  that  the  element  of  solid  angle  is  sin0ed0B 
(the  <{>B  integration  has  already  been  performed),  it  is  found  that 

if  =  ko  1 2 ( Kof  -  KjJfl  +  kgsf(Kof  -  Kj)]'! 

+  ( K0f  -  K3)(l  +  k0Sf(K0f  -  K3)]-1j/3  (38) 

With  either  of  these  expressions  for  in  (34)  the  quasi-static 

ts 

dielectric  constant  may  be  found.  A  numerical  solution  is  necessary. 

V.  SCATTERING  TERMS 

The  scattering  contribution  E  to  the  effective  dielectric  con- 

% 

stant  K  in  (14)  will  now  be  considered.  Scattering  within  a  grain  is 
neglected  because  of  the  very  small  diameters  and  close  spacing  of  the  brine 
cells  within  the  qrain  compared  to  the  wavelengths  of  interest.  Thus,  the 
sea  ice  is  treated  as  a  two-component  material  consisting  of  grains  with 
dielectric  constant  (32)  and  air  bubbles.  The  scattering  correction  depends 
only  on  scatterina  from  the  air  bubbles  and  scatterina  from  the  grains  due  to 
the  fact  that  the  dielectric  tensors  describing  the  individual  grains  have 
randomly  oriented  principal  axes  corresponding  to  randomly  oriented  grain 
c  axes. 

There  are  three  factors  to  be  considered  in  evaluating  (15).  These 
are  the  specifications  of  C^jj^  »  ,  and  the  evaluation  of  the  principal 

value  integral.  First  consider  the  fourth  rank  correlation  tensor  (16). 
Following  the  method  discussed  in  (15],  £  is  written  as 
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Ur)  =  4air  H(r)  +  £9  [1  -  H( r )  1 

*V 

=  lair  H(r)  +  L  11  -  H(  r )  ]  +  {£q  -4q)[l  -  H(r)  ]  (39) 

w  ~  as  ~  ~  ~  ^  ~ 


where  H(r)  is  a  random  function  equal  to  one  if  r  is  in  an  air  bubble  and 
zero  if  r  is  in  a  grain.  When  forming  C^jj^  the  first  two  terms  in  the 
second  line  of  (39)  are  evaluated  exactly  as  in  (15).  The  third  term  does 
not  lead  to  cross  products  with  the  first  two  terms  when  expected  values  are 
taken  so  that 


cijk£(£‘E’ ) 


In  (40) 


Bijk^r2“Il> 


t4air ij  4air k2  +  4gij  4qk£  “  4air ij  4gk£ 

~  fgij  4airkglRa<£  “  £'  >  +  Bijk2<£2  ~  £l  >  <40> 

4g  i  j d  )4g  k£  ^ ^ )  -  4g  ij  C  g  kf  111  -  H(rj)][l  -  H(£2))^ 

(41  ) 


with  the  arguments  1  and  2  in  |q  referrina  to  the  Euler  anqles  and  Ra  is 

ss 

a  scalar  function  describing  the  air  bubble  statistics  (15].  It  will  be 
assumed  that  the  air  bubbles  have  an  exponential  correlation  function  so  that 
[15] 

Ra(£2~£l>  =  vair(1  ~  vair>  exP! " I £2"£l| a )  (42) 

where  the  correlation  lenqth  £a  is  related  to  the  mean  bubble  diameter  da 
by 

«a  =  f  (1  ”  vair ^  da  (43) 

The  term  in  (40)  is  considerably  more  complex  and  much  data  on  ice 

remains  to  be  analyzed  to  determine  its  precise  behavior.  However,  it  is 
obvious  that  Bjj^  approaches  zero  as  |r2“£l|  becomes  laroe.  The  ad-hoc 
assumption  will  be  made  that 
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Bijk$(l2~£l)  *  Bijk£<°>  exPf-|£2'£l|/^qI  (44> 

From  (40),  is  evaluated  as 

Bijk£*0^  =  ij  £g  k  £  ~  4g  ij  £g  k0lvgrain  (45) 

According  to  the  overlap  picture  discussed  in  [15]  it  is  expected  that  the 
correlation  length  will  be  somewhat  less  than  the  mean  grain  diameter 

dg  .  It  will  be  assumed  to  be  qiven  by 

tq  =  f  dg  (46) 

The  second  term  in  (45)  has  already  been  evaluated  in  (36)-(37). 
Before  writing  equations  for  the  first  term,  it  is  useful  to  note  the  sym¬ 
metries  of  Bjjkt'1  .  Bijk(?(°)  is  symmetric  under  the  interchange  of  i  and 

j  ,  k  and  £  ,  and  the  pair  (ij)«»—  (k£)  as  can  be  seen  by  examination  of  (36) 

Further,  the  basic  assumption  of  invariance  under  orthogonal  transformations 
leaving  the  z  axis  fixed  holds-  All  of  these  symmetries  combined  reduce 
the  total  number  of  independent  components  to  only  five-  In  fact,  it  is 
found  that 

Bijk2<0)  =  al  &ij  bk£  +  a2 1  ^ik  6j  g  +  &i£  6  jk  1  +  bl^ij  &k3  ^3 

+  k£  ^i3  ^  j  3  ^  +  b2^ik  ^  j3  &k3  +  6 il  6j3  &k3  +  6jk  &i3  bi3 
+  &i3  6 k 3 1  +  <=i  &i3  6 j 3  6*3  &£ 3  <47) 

where  the  6’s  are  the  Kronecker  deltas  and  the  constants  are 


al 

=  B1122(0> 

(48a) 

a2 

=  Bj  21  2  ( 0  > 

(  48b) 

bl 

=  B1133(0)  -  Bj  2  22 ( 0  ) 

(  48c  ' 
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b3  =  B3333(0)  “  B1212^®^  ( 48d 

c3  *=  ®3333^0)  -  a3  -  2a2  -  2b3  -  4b2  (48e 

A  routine  but  lengthy  computation  of  the  average  of  the  product  in  (45)  may 
be  performed.  The  average  over  <f>B  is  simply  an  integral  and,  except  for 
frazil,  the  0B  average  is  treated  as  in  Section  IV  by  interpreting  0B  to 
be  a  parameter  representing  a  mean  brine  cell  axis  angle.  For  simplicity, 
this  angle  is  identified  with  the  parameter  0B  appearing  in  Section  IV. 
With  these  approximations  the  result  is 

b1122(0)  =  ,*4  fs302  cos20b  sin2eB/{olY)]2  _  s3p2  cos26b  sin2eB/ 


K 

( 48d ) 

m 

>  V 

(48e) 

*»  <i 

*  V  ^ 

may 

K 

-  |[psin20B(o3  +  S3p)/(o,y) ]2  J  vgrain 
b1212(0)  =  ^  Jp  sin20B(O3  +  S3  J  2  vgrain 

b1313(0)  =  \  [ p  cos6B  sin0B/ y] 2  vgrain 
B1133(0)  =  B3333(0)  *  0 


(49a) 

(49b) 

(49c) 

( 49d ) 


where  the  quantities  that  appear  have  been  defined  in  equations  (34)  —  38). 

As  in  the  case  of  the  quasi-static  dielectric  constant,  frazil  may 
be  treated  more  precisely.  In  fact,  the  spherical  symmetry  implies  that  b3 
=  b2  =  C)  =  0  in  (47).  Contracting  on  the  two  inner  and  two  outer  indices 
in  (47)  shows  that 

Bi  j  j i ( 0  >  =  3al  +  12a2  (50) 

But  recalling  that  K0  ,  S  ,  and  £g  are  scalar  multiples  of  the  unit  dyad 

ss  ss  ~ 

for  frazil  allows  Bijji(0)  to  be  evaluated  using  (36)  as 

Bijji^*  =  ^o  vgrain  j  2  ^  K0f“  Kl)^tl  +  kg  Bf^Kof~  K1  ^  1  2 

+  <Kof-  K3>2n  +  kg  Sf(Kof  -  K  3  )  )  “ 2  j  -  3j{2  varain  (51) 
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Further,  a2  *n  (47b)  rosy  be  evaluated  exactly  by  integrating  (49b)  after 
multiplying  by  —  sin8g  «  The  result  is 

a2  -  ip/{o-i^3)  l2  vgrain/15  (52) 

where  and  ©3  are  to  be  specialized  to  the  frazil  parameters. 

The  next  problem  in  evaluating  (15)  is  the  Green's  function  which 

satisfies  (9)  when  K0  has  the  form  (17).  It  is  not  difficult  to  find  the 

Fourier  transform  of  F  [17].  An  evaluation  of  the  inverse  transform  yields 

as 

r  ( r )  =  (4  it  kg  K01  Kq32)  ~17V(e1/rk)  +  Fj  (53) 

as  ~  '  555 

where  the  nonvanishing  components  of  are 


<ri»n 

= 

eD/  ( 4-rrr )  -ojy2  -  x2  +  x2 

l el  -  eol 

(54a) 

(rih2 

S 

e0/(4irr)  -  ojx2  -  y2  +  y2 

lei  -  eD] 

(54b) 

(ri)i2 

= 

<r 1  )2i  =  x y 0 [2  -p^ 

lei 

~  e0) 

(54c) 

(ri)33 

= 

1/2 

e1/(4rK03  rK ) 

( 54d ) 

p 

= 

[x2  +  y2 ] 1/2 

(55a) 

rK 

= 

[p2/K01  +  z2/K02 ] 1/2 

(  55b) 

O 

= 

i/(*T  kQ  Kj{2  P«) 

(55c) 

eo 

= 

1/2 

exp[ i  kQ  K01  r] 

(  55d  ) 

el 

= 

exp(i  ko(K01K03)1/2rK] 

(  55e ) 

Finally,  turninq  to  the  integral  in  (15),  note  that  (53)  shows 

that  r*  is  a  symmetric  tensor.  This  combined  with  (40),  (47),  and  (4R) 
% 

shows  that  the  tensor  =  C^jk£Fjk  is  symmetric.  Its  components  are 


17 
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lc1122  +  c1212 i  +  Ic1313  "  c1212lr33  +  c1212  Tr  £' 


r C1 1 22  +  c1212lr22  +  f°1313  ~  c1212)r33  +  c1212  Tr  £' 


(56a) 


(56b) 


t c3333  "  c1313lr1313  +  c1313  Tr  £' 


(56c) 


P21  =  f c1122  +  c1212Jri2 


( 56<i  ) 


=  P31  =  I c1133  +  c1313lr{3 


( 56e ) 


p23  =  p32  =  fc1133  +  c1313 JT23 


(56f  ) 


where  the  symbol  Tr  means  trace. 

Thus,  in  view  of  the  assumed  exponential  correlation  functions  (42) 

and  (44),  an  evaluation  of  inteqrals  involvinq  components  of  r*  and  exponen¬ 
ts 

tial  functions  are  necessary  to  evaluate  E  .  We  note  immediately  that  terms 
involvinq  will  vanish  because  of  the  factor  y  in  (54c)  and  the  fact 

that  k  has  no  y  component  (see  (18)).  Similarly,  the  V  V  operation  in 
(53)  will  lead  to  E32  =  E23  =  0  .  These  observations  merely  confirm  the 
form  (19)  already  deduced  for  K  .  For  the  remainino  components,  a  prelimi- 
nary  inteqration  by  parts  for  those  terms  arisinq  from  VV(e^/rK)  in  (53) 
is  useful  to  allow  an  easy  evaluation  of  the  principal  value  required  in  (15). 
Letting  C(r)  be  either  of  the  exponentials  in  (42)  or  (44)  and  recallinq 
(11),  it  is  found  that 

Pvjexpl -i  k- r  ]C( r )?V( ej/rK )dV  =  j exp[ -i  k -r  ]C ( r )VV ( e3/rK )dV  -  _S 


=  |  V  jexp( -i  k -r  ]C(r  )j  V(ej/rK)dV  -  ^S 


because  the  inteqral  over  the  surface  at  infinity  vanishes.  Since  the 
singularity  of  V(e^/r^)  is  sufficiently  mild,  no  particular  precautions 
need  be  taken  near  the  origin  in  evaluatinq  the  resultant  inteoral. 
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Aside  from  easy  inteqrals  involvinq  the  terms  e0/(4Trr)  in  (54a) 
and  (54b),  it  has  not  been  found  possible  to  evaluate  the  remaining  integral 
for  E  in  closed  form.  However,  if  spherical  coordinates  r  ,  6  e  4>  are 
introduced,  the  4>  integration  leads  to  Bessel  functions.  The  inteqrals  of 
these  with  exponentials  in  r  over  the  range  from  zero  to  infinity  are  tabu 
lated  inteqrals.  Thus,  only  integrals  over  the  range  0  to  tt  in  0 
remain  to  be  evaluated  numerically.  Because  of  their  lenqth,  the  formulas 
will  not  be  written  here  although  the  equations  already  discussed  make  their 
derivation  routine.  They  do  not  contain  singularities  and  may  be  evaluated 
using  Gauss  quadrature  formulas.  The  most  difficult  numerical  situation 
arises  at  low  frequencies  for  columnar  ice  of  high  salinity  where  Kqj  is 
siqnif icantly  different  from  Kq3  .  It  has  been  found  that  a  10-point 
quadrature  will  result  in  errors  much  less  than  one  part  in  a  thousand  for 
the  components  of  K  at  frequencies  as  low  as  0.1  GHz. 
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VI.  COMPARISONS  WITH  EXPERIMENTAL  DATA 

Most  reported  experiments  on  the  dielectric  constant  of  sea  ice 
have  treated  a  situation  corresponding  to  a  vertically  propagating  wave  inci¬ 
dent  on  a  slab  of  ice  oriented  in  its  natural  growth  position.  In  the  nota¬ 
tion  of  this  paper,  kx  =  0  in  (18).  Thus  the  K22  component,  which  is 
equal  to  the  K33  component  because  kx  =  0  ,  of  the  tensor  K  was  measured 
In  a  few  cases  [2],  (18]  additional  measurements  with  the  polarization  of  the 
electromagnetic  wave  parallel  to  what  is  called  the  z  axis  in  this  paper 
were  made.  Thus  kz  =  0  and  the  K33  component  was  measured.  In  [2], 
further  measurements  were  performed  on  ice  samples  cut  at  acute  angles  to  the 

growth  direction.  Some  complicated  function  of  the  components  of  K  was 

a 

measured  because  both  ordinary  and  extraordinary  waves  would  be  generated  in 
the  sample.  Unfortunately,  not  enough  details  concerning  the  measurements 
were  presented  in  these  cases  to  allow  a  detailed  analysis. 

Experimental  data  from  (1]  on  two  samples  of  first-year  ice  denoted 
by  FY  1-4  and  FY  1-9  are  considered  first.  Fiaures  1  and  2  show  comparisons 
of  computations  based  on  the  theory  of  this  paper  and  measurements  of  the 
real  and  imaginary  parts  of  K22  at  two  distinct  temperatures.  The  computa¬ 
tions  use  the  salinity  S  and  density  p  given  in  [1]  together  with  an 
assumed  mean  brine  cell  angle  83  =  24°  ,  correlation  length  ratio  @z/fp  = 
200  ,  a  mean  air  bubble  diameter  da  =  0.49  mm  ,  a  mean  ice  grain  diameter 
dg  =  1  cm  ,  and  the  expression  (30)  for  Sgj  ,  which  is  based  on  a  cylindrical 
correlation  function.  If  (27)  is  used  for  Sqi  ,  the  computed  real  parts  are 
found  to  be  slightly  larger  near  0.1  GHz  and  almost  unchanged  at  1  GHz  but 
the  computed  imaginary  parts  are  significantly  larger  (e.g,,  1.7  versus  1.48 
for  FY  1-9  at  t  =  -5 . 2 °C )  at  0.1  GHz.  Even  at  1  GHz,  the  use  of  (30) 
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instead  of  (27)  results  in  a  noticeable  improvement  of  the  fit  for  Im  K22  . 
It  is  concluded  that  the  low  frequency  data  does  not  support  the  assumption 
of  an  elliptical  equi-correlation  surface  for  the  brine  cells.  A  comparison 
of  the  results  shown  in  Figures  1  and  2  with  the  theories  discussed  in  [1] 
and  [5]  show  a  considerably  improved  fit  to  the  data.  This  is  achieved  with 
a  more  realistic  brine  cell  length  to  diameter  ratio.  In  fact,  a  ratio 
£z/#p  =  20  such  as  used  in  [1]  and  (5]  yields  a  poor  fit  to  the  experimental 
data . 

Although  K33  was  not  measured  for  these  two  samples,  an  examina¬ 
tion  of  computed  values  of  K33  with  k2  =  0  is  informative.  It  is  found 
that  K33  is  considerably  larger  than  K22  in  the  frequency  range  of 
Figures  1  and  2.  For  example,  K33  =  29. 25+41. 74i  at  0.1  GHz  and  5.64+3. 37i 
at  10  GHz  for  the  FY  1-9  sample  at  -5.2°C.  Even  for  the  FY  1-4  sample  at 
-1 4 . 8 °C ,  K33  =  10.54+9.29i  at  0.1  GHz  and  3.48+0.51i  at  10  GHz.  At  still 
higher  frequencies.  Re  K33(kx  =  0)  approaches  Re  ^23(^7.  =  0)  not  only  for 
these  cases  but  for  all  sea  ice  according  to  computations.  For  ice  with  a 
salinity  greater  than  a  few  tenths  ppt  (parts  per  thousand),  however,  the 
imaginary  parts  of  the  components  of  K  differ  measurably. 

These  computations  are  verified  by  a  number  of  observations.  Thus 
Sackinger  and  Byrd  [18]  have  measured  the  dielectric  constant  of  ice  samples 
in  the  26-40  GHz  range.  As  an  example,  we  consider  their  measurements  on  ice 
with  a  salinity  S  =  7.2  ppt  in  the  temperature  range  -16.5  to  -32°C.  Essen¬ 
tially  identical  results  were  found  for  the  real  part  of  the  dielectric  con¬ 
stant  whether  the  polarization  of  the  electric  field  was  perpendicular  or 
parallel  to  the  qrowth  direction  (K22  with  kx  =  0  and  K33  with  kz  =  0 
in  the  present  notation).  The  ice  density  was  not  stated  in  (181.  However, 
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if  p  =  0.83  g/cm3  is  assumed,  the  computed  Re  K^2  values  at  26  GHz  range 
from  3.12  at  -16.5°C  to  2.90  at  -32°C  in  agreement  with  the  reported  values. 
Further  the  computed  Re  K33  is  larger  than  these  values  by  only  0.05  and 
0.008  at  -16.5  and  -32°C  respectively.  Also,  the  computations  show  only  an 
insignificant  frequency  variation  of  Re  K  in  the  26-40  GHz  range.  Again, 

S3 

this  is  verified  by  the  data  in  [18].  The  situation  is  different  for  Im  K 
and  is  summarized  in  Table  1.  Unfortunately  there  was  considerable  scatter 
in  the  reported  measurements  so  that  only  mean  values  and  an  indication  of 
the  scatter  in  the  data  is  given.  Computed  values  in  the  range  26-40  GHz 
are  also  shown  in  the  table  where  the  first  number  corresponds  to  26  GHz  and 
the  second  to  40  GHz.  In  addition  to  the  above  stated  parameters,  the  param¬ 
eters  used  in  these  computations  are  0B  =  24°  ,  £z/£p  =  200  ,  da  =  0.5  mm.  , 
and  da  =  1  cm  .  It  is  seen  that  the  computations  generally  track  the 
observations . 

Other  relevant  experiments  are  discussed  in  [2].  Measurements  on 
old  ice  with  S  s;  0  ppt  at  a  temperature  of  -13°C  were  made  at  frequencies  of 
3.79,  10,  and  38.46  GHz.  Because  of  the  low  salinity  the  measurements  of 
Re  K22  and  Re  K33  (corresponding  to  4>  =  90°  and  4>  =  0°  respectively  in 
the  notation  of  [2])  showed  essentially  no  anisotropy.  Values  were  obtained 
in  the  ranqe  2.5  to  2.7.  Althouqh  no  density  information  was  qiven  in  [2], 
computations  assuming  p  =  0.7  a/cm3  ,  which  is  reasonable  for  old  ice, 
yielded  values  of  Re  K22  ranqinq  from  2.51  to  2.59  between  3.79  and  38.46 
GHz  for  S  =  0.05  ppt.  The  larqest  difference  between  Re  K33  and  Re  K22 
was  0.005.  These  computations  also  assumed  0B  =  24°  ,  t^/^p  =  200  ,  dg  = 

1  cm  ,  and  da  =  1.2  mm  .  The  assumed  mean  air  bubble  diameter  is  quite 
reasonable  for  old  ice  and  yielded  an  attenuation 
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a  =  20  kQ  log10e  In  K1/2  (58) 

(where  K  may  be  either  K22  or  K33)  of  2.01  dB/cm  when  kx  =  0  and  2.04 
dB/cm  when  kz  =  0  for  a  frequency  of  38.46  GHz.  This  is  in  agreement  with 
the  measured  value  of  2  dB/cm  for  each  of  the  orientations.  At  38  GHz,  Im  K22 
is  over  33  times  greater  than  the  quasi-static  Im  Kq22  so  that  the  scatter¬ 
ing  contribution  to  K  is  responsible  for  essentially  the  entire  measured 
attenuation  of  this  sample  of  ice.  The  attenuation  at  3.79  and  10  GHz  was 
too  low  to  be  measured.  Computations  show  a  small  anisotropy  at  10  GHz.  The 
computed  result  for  the  K33  component  is  0.05  dB/cm  at  10  GHz  which  would 
appear  to  be  difficult  to  measure  based  on  the  indicated  system  sensitivity  [2] 
Measurements  on  two  samples  of  winter  sea  ice  with  salinities  of 
6.55  and  6.2  ppt  at  temperatures  of  -13  and  -12°C  respectively  were  also  dis¬ 
cussed  in  [2].  These  samples  showed  considerably  different  behavior  than  the 
pack  ice.  Again,  the  ice  densities  were  not  given.  Based  on  the  rather  larae 
reported  values  of  Re  K  ,  the  assumption  was  made  that  p  =  0.9  g/cm^  with 
the  other  parameters  required  by  the  theory  taken  as  0B  =  24°  ,  =  200  , 

da  =  0.5  mm  ,  and  dq  =  1  cm.  Comparisons  with  the  measurements  are  shown  in 
Tables  2  and  3  for  the  two  ice  salinities.  It  is  seen  that  the  measured  real 
parts  are  higher  than  computed  values  at  3.79  and  10  GHz.  However,  the  meas¬ 
ured  difference  Re  K33~Re  K22  and  the  calculated  difference  Re  1(33-80 
agree  closely  except  for  10  GHz  in  Table  3  where  no  anisotropy  is  shown  in 
the  measurement.  The  hiqh  measured  values  do  not  appear  to  be  consistent  with 
the  data  in  Figure  1  at  4  GHz  and  -14.5°C  where,  for  a  much  higher  salinity, 
a  lower  value  of  Re  K22  was  obtained.  The  calculated  attenuation  values 
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also  show  some  points  of  agreement  with  the  measurements  although  there  are 
dif ferences.  There  are  too  many  unknowns  to  speculate  as  to  the  origin  of 
these  differences. 

Another  source  of  data,  which  is  interesting  because  frazil  is 
identified  explicitly,  will  be  examined.  In  [3],  measurements  at  a  frequency 
of  10  GHz  were  made  on  two  samples  of  frazil  as  a  function  of  temperature. 
Figures  3  and  4  show  computed  and  measured  results  for  Re  K  and  Im  K  . 
According  to  the  discussion  in  Sections  IV  and  V,  the  propagation  constant  is 
independent  of  direction  and  polarization  but  the  experiment  did  not  attempt 
to  verify  this  property.  The  computations  assume  £z/£p  -  200,  da  =  1  mm,  dg 
=  4  mm  in  addition  to  the  densitites  and  salinities  stated  in  [3].  Notice  that, 
at  some  temperatures,  two  measured  values  are  shown.  These  correspond  to  meas¬ 
urements  made  during  increasing  and  decreasing  temperature  cycles  and  provide 
some  indication  of  the  random  experimental  errors.  Figure  4  shows  that  rather 
close  agreement  with  the  measured  Im K  values  is  achieved.  Figure  3,  however, 
indicates  a  consistent  difference  of  about  0.2  between  the  measured  and  com¬ 
puted  values  of  Re  K  over  most  of  the  temperature  range.  This  difference  is 
somewhat  surprising  -  especially  at  temperatures  below  -20°C  where  the  brine 
volume  in  ice  of  this  salinity  becomes  negligible.  The  measured  Re K  values 
approximate  those  of  pure  solid  ice,  not  porous  ice  with  a  density  near  0.84 
g/cm^  which  has  a  dielectric  constant  whose  real  part  is  near  2.87.  On  the 
other  hand,  the  hydrated  salts  in  the  ice  may  have  some  unexpectedly  larqe 
effects  which  are  not  treated  as  being  substantially  different  from  ice  in 
the  theory. 

Data  at  10  GHz  on  columnar  and  multiyear  ice  are  also  discussed  in 
(3].  The  two  columnar  ice  cases  are  puzzlinq  and  not  in  accord  with  physical 
expectations.  The  ice  with  the  higher  density  and  higher  salinity  had  a 
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measured  Re  K22  much  lower  than  that  of  the  lower  density,  lower  salinity 
ice.  Even  its  Im  K22  value  was  lower  for  temperatures  above  -30°C.  Both 
the  higher  salinity  and  higher  density  should  have  had  the  opposite  effect. 

For  the  two  multiyear  ice  cases  in  [3]  the  computed  values  of  Re  K22  are 
essentially  the  same  and  higher  than  the  measured  values  by  about  0.2.  Lower¬ 
ing  the  assumed  ice  density  to  0.73  g/cm^  instead  of  the  specified  0.77  g/cm^ 
for  the  S  =  0.61  ppt  case  brings  the  calculated  values  into  agreement  with 
the  measurements  at  all  temperatures.  The  calculated  imaginary  parts,  however 
are  low  by  a  factor  of  about  0.5.  This  may  not  be  a  problem  because  the 
scatter  in  the  measured  values  indicate  that  high  accuracy  was  not  achieved 
in  these  measurements  for  such  low  Im  K22  values. 

VII.  CONCLUSIONS  AND  RECOMMENDATIONS 

A  theory  of  the  dielectric  constant  of  sea  ice  has  been  developed 
which  is  able  to  account  quantitatively  for  a  large  amount  of  published  data 
in  the  frequency  range  0.1  to  40  GHz. 

However,  some  disagreements  were  noted.  These  may  be  due  to  a  lack 
of  quality  in  the  data  for  there  are  apparent  problems  among  some  of  the 
measurements  themselves.  The  lack  of  really  precise  measurements  shows,  for 
example,  in  the  data  discussed  with  respect  to  Figures  1  and  2  where,  even 
though  the  theory  clearly  reproduces  the  main  observed  characteristics ,  the 
experimental  points  do  exhibit  what  must  be  assumed  to  be  some  significant 
random  errors  (the  assumption  that  these  correspond  to  real  changes  appears 
to  be  quite  unreasonable  and  no  known  physical  process  could  account  for  them). 
The  large  scatter  in  the  experimentally  determined  imaainary  part  of  the  die¬ 
lectric  constant  data  discussed  in  connection  with  Table  1  was  noted  pre¬ 
viously.  This  precluded  a  more  quantitative  check  of  the  theory,  which 
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generally  falls  within  the  reported  range  at  each  temperature,  as  a  function 
of  frequency  for  the  particular  data  source  in  question.  Bias  type  errors 
are  also  of  some  significance.  For  example  the  Re  K22  measurements  dis¬ 
cussed  in  connection  with  Tables  2  and  3  at  3.79  and  10  GHz  show  a  definite 
bias  with  respect  to  experimental  data  on  approximately  similar  ice  (FY  1-4) 
shown  in  Figure  1.  The  correction  of  this  bias  at  3.79  and  10  GHz,  if  it  is 
real  and  the  same  for  Re  K33  as  for  Re  K22  »  would  have  resulted  in  sig¬ 
nificant  improvement  in  agreement  between  the  calculated  and  measured  values 
of  Re  K33  as  well  as  Re  K22  - 

These  observations  clearly  show  that  dielectric  measurement  tech¬ 
niques  for  this  difficult  to  measure  substance  must  be  developed  further. 
Complete  measurements  require  data  not  only  for  the  case  where  the  electric 
field  polarization  is  perpendicular  to  the  ice  growth  direction,  but  also 
parallel  to  it.  As  this  report  has  shown,  there  is  ample  evidence  that  most 
forms  of  sea  ice  have  anisotropic  dielectric  properties.  Measurements  at 
several  frequencies  and  temperatures  on  a  given  sea  ice  sample  should  be 
mandatory.  This  procedure  would  greatly  facilitate  the  comparison  and  check¬ 
ing  of  data  from  diverse  sources  since  the  large  number  of  factors  which 
influence  the  dielectric  properties  of  sea  ice  would  otherwise  make 
intercomparison  very  difficult. 

In  addition  to  the  dielectric  constant  measurements  themselves, 
more  care  should  also  be  exercised  in  reporting  on  such  obvious  physical 
properties  such  as  density,  salinity,  and  temperature.  These  play  a  major 
role  in  accounting  for  the  dielectric  properties  even  though,  as  seen  in  the 
discussion  of  the  existing  experimental  data,  they  are  not  always  oiven.  A 
more  subtle  area  is  the  specification  of  micro-structure  statistics  on  such 
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features  as  air  bubble  size  and  brine  cell  length,  dimension  and  orientation. 
These  have  not  customarily  been  measured  in  connection  with  dielectric 
studies  in  the  past  but  do  seem  to  play  a  role.  In  fact,  without  these 
latter  parameters,  there  is  little  hope  of  refining  the  theory  since  these 
have  been  shown  to  influence  the  dielectric  properties.  The  lack  of  such 
published  data  required  a  number  of  ad-hoc  assumptions  in  the  present  study. 
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Imaginary  part  of  X  for  frazil  at  10  GHz 
as  a  function  of  temperature. 
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TABLE  1  COMPARISON  OF  CALCULATIONS  OF  THE  IMAGINARY  PARTS  OF  THE 
COMPONENTS  OF  K  WITH  MEASUREMENTS  ON  SEA  ICE  BETWEEN 
26  AND  40  GHZ  (S  =  7.2  ppt) 


1 

Temperature 

(°C) 

Im  K 

22 

Im  K33 

Calculated 

Observed 

Calculated 

Observed 

-32.0 

0.0168  to  0.0244 

0.011  +0.002 

0.0397  to  0.0388 

0.031  +0.003 

-25.0 

0.0348  to  0.0409 

0.045  +0.015 

7 

0.110  to  0.087,8 

0.17  +0.06 

-21.5 

0.0626  to  0.0668 

0.054  +0.025 

0.225  to  0.168 

0.19  +0.06 

-16.5 

0.0780  to  0.084 

0.084  +0.036 

0.320  to  0.237 

0.22  +0.06 

TABLE  2  COMPARISON  OF  WINTER  ICE  (S  -  6.55  ppt,  t  =  -13°C) 
DATA  IN  [21  WITH  THEORY 


Frequency 
(GHz  ) 

|  Measured 

Calculated 

Re  K33 

Re  K22 

a(dB/cm)  j 

Re  K33 

Re  K22 

Or  ( dB/cm ) 

m 

kx  =  0 

rasa 

kx  =  0 

3.79 

4.25 

3.55 

0.90 

0.25 

4.06 

3.34 

1.79 

0.20 

10.0 

3.85 

3.45 

2.8 

1.0 

3.63 

3.29 

3.54 

0.43 

38.46 

3.13 

3.13 

5.6 

2.2 

3.27 

3.23 

5.68 

1.64 

_ 

Measured 


Calculated 


Frequency 
(GHz  ) 

O'  ( dB/cm ) 

0 ( dB/cm ) 

Re  K33 

Re  K22 

am 

EBB 

Re  K33 

Re  K22 

kz  =  0 

BBS 

3.79 

4.6 

3.9 

1.1 

0.2 

4.01 

3.32 

1.71 

0.19 

o 

o 

H 

3.8 

3.8 

3.0 

1.2 

3.60 

3.29 

3.37 

0.42 

38.46 

3.1 

3.2 

5.2 

2.8 

3.26 

3.23 

5.41 

. 

1.56 

